Molecular Genetics and Genomics

BackgroundConventional models of human ribosomal DNA (rDNA) array organization have historically depended on transcription-centric boundaries, partitioning the unit into a [~]13 kb rDNA transcription region and a monolithic [~]31 kb intergenic spacer (IGS). While our previous identification of Duplication Segment Units (DSUs) mapped these arrays based on an intuitive analysis of the microsatellite density landscape of the complete reference human genome, our present deep mining of this landscape has revealed a more accurate rDNA Gene Unit Pattern. Methods & ResultsIn this study, we conducted a deep mining analysis of our previously established microsatellite density landscape of the T2T-CHM13 assembly, focusing specifically on nucleolar organizing regions (NORs). We suggest a more accurate rDNA Gene Unit Pattern containing a (CTTT)n microsatellite aggregation ahead of the rDNA gene and a (CT)n microsatellite aggregation behind the gene, rather than a pattern featuring an IGS region inserted between two rDNA genes. ConclusionsA correct rDNA gene pattern of the human genome probably includes a (CTTT)n microsatellite aggregation ahead of the gene and a (CT)n microsatellite aggregation behind it, which possibly constitute cis- and trans-regulating regions; the (CTTT)n and (CT)n microsatellite aggregations may provide two different local stable DNA structures for regulatory protein binding.

BackgroundPolycystic ovary syndrome (PCOS) is a prevalent endocrine disorder and a leading cause of female infertility, with complex genetic, metabolic, and hormonal etiologies. Long non-coding RNAs (lncRNAs) have emerged as important regulators of diverse biological processes, yet their roles in PCOS remain underexplored. Here, we identified and characterized PCOS differentially expressed gene-associated lncRNAs (PDEGAL) with an integrative approach combining expression data, genetic association, and evolutionary analysis. MethodsThirty-three PCOS-associated protein-coding genes were obtained from our prior study, and all their nearby and overlapping lncRNAs were annotated. These candidates were analyzed using UCSC Genome Browser-mapped annotations and datasets, including NCBI RefSeq, GENCODE, GTEx, GWAS SNPs, and conservation, as well as the FANTOM5 cap analysis of gene expression (CAGE) promoter data, to assess their expression, regulatory potential, genetic variant overlaps, and evolutionary conservation. ResultsTwenty-three PDEGALs (18 antisense to, and 5 sharing bidirectional promoters with, known PCOS-associated protein-coding genes) were identified. 17 PDEGALs contained GWAS SNPs with statistically significant disease associations, 9 of which were associated with PCOS-related traits. 5 PDEGALs demonstrated expression in the KGN granulosa cell model of PCOS. Key gene structure element (KGSE) analysis revealed that most PDEGALs are primate-specific. Integrating four criteria--GTEx expression, GWAS SNPs, FANTOM promoterome, and KGSE conservation--highlighted HELLPAR as the only lncRNA fulfilling all four, while five others--PGR-AS1, MTOR-AS1, ENSG00000265179, ENSG00000256218, and LOC105377276--fulfilled three of the four criteria. ConclusionsWe have systematically identified candidate PCOS regulatory lncRNAs with convergent genetic, expression, and evolutionary evidence. These results provide a framework for functional validation and highlight lncRNAs as potential biomarkers and therapeutic targets in PCOS that function by regulating their nearby and overlapping protein-coding genes.

Agriculture has modified the soil structure due to the influence of external factors and processes that affect microbial biodiversity. Metagenomics is a fundamental tool for the study of soil microbial diversity because it provides information about the ecosystem diversity, including both the microorganisms that cannot be isolated in culture media and those that are no longer viable in the analyzed sample. In this work, six soil samples obtained from agroecosystems of central and northern Argentina were subjected to a preliminary 16S metagenomic analysis. Copiotrophic bacteria (Proteobacteria and Actinobacteria) were dominant and one of the samples had a dominance of an oligotrophic Phylum (Acidobacteria). Our findings support previous evidence from traditionally managed agroecosystems and provide new insights into the diversity of soil microbiomes in Argentine regions outside the Pampas. Finally, we analyzed the most common genera with relevant species to agronomy, both beneficial and pathogenic, and their abundance and diversity in the sequenced samples.

Wellbeing is commonly defined as the combination of feeling good and functioning well and typically conceptualized as two related but distinct components. Hedonic wellbeing emphasizes pleasure, happiness, and life satisfaction, while eudaimonic wellbeing focuses on meaning, personal growth, flourishing, and the realization of ones potential. The Mental Health Continuum-Short Form was developed as a comprehensive measure of wellbeing and includes three subscales assessing emotional, social, and psychological wellbeing. Although the Mental Health Continuum total score is often interpreted as an indicator of overall wellbeing, the underlying genetic structure of its three subscales and its genetic overlap with other commonly used wellbeing measures remains unclear. Using data from 5,212 individuals from the Netherlands Twin Register (72% female, mean age 36.4), we fitted multivariate twin models to examine the genetic architecture of the Mental Health Continuum and its associations with other wellbeing measures (quality of life, life satisfaction, subjective happiness, and flourishing). Results indicate that, at the genetic level, the Mental Health Continuum is best explained by its three distinct subscales rather than by a latent factor. When considering the Mental Health Continuum together with the other wellbeing measures, we found moderate to high genetic correlations (r = 0.52 - 0.83), indicating substantial overlap in the genetics underlying the wellbeing constructs. However, we did not find evidence for a single common genetic factor underlying all constructs. These findings highlight the multidimensional structure of wellbeing, but the moderate to high genetic correlations across measures suggest that it is important to align the level of measurement (phenotypic vs genetic) with the research question.

Psychosocial stress is a significant risk factor for mental and physical illness, and emerging evidence suggests that altered oral microRNAs (miRNAs) and microbiome may act as biomarkers or mediators of stress responses. This study investigated stress-associated molecular changes in saliva from 113 male police officers. Based on repeated administrations of the Karasek Demand/Control and Effort/Reward Imbalance questionnaires, subjects were stratified by perceived stress response (SR) to homogeneous occupational stressors into low, intermediate, or high responders. Salivary miRNA profiles were analyzed using small RNA sequencing, and microbiome composition was assessed through shotgun metagenomics. Eighteen miRNAs were significantly differentially expressed between high- and low-SR groups, with four miRNAs with increasing (miR-10400-5p, miR-1290, miR-6074-5p, and miR-9902) and fourteen with decreasing (including miR-21-5p and mirR-142-3p) levels in the high SR group (adj.p<0.05). The identified salivary miRNAs showed a progressive alteration from low- to high-SR groups. Functional enrichment analysis indicated that dysregulated miRNA targets are involved in apoptosis, cellular stress responses, and metabolic regulation. Distinct salivary microbial communities were also observed across SR groups. Several taxa displayed progressive abundance shifts, with Prevotella baroniae and Schaalia odontolytica increasing and Actinomyces naeslundii and Capnocytophaga ochracea decreasing in the high SR group. Functional predictions revealed, in this group, a significant enrichment of inositol degradation pathways, paralleled by a reduction in bacteria involved in L-tryptophan and thiamine biosynthesis. These findings suggest that salivary miRNAs and microbiota profiles may serve as non-invasive biomarkers of psychosocial stress and provide insight into molecular mechanisms linking chronic stress to physiological and behavioral outcomes.

Proteins and non-coding RNAs are functional products of the genome that are central for crucial cellular processes. With recent technological advances, researchers can sequence genomes in the thousands and probe numerous genomic activities of many species and conditions. Such studies have identified thousands of potential proteins, RNAs and associated activities. However there are conflicting interpretations of the results and therefore which regions of the genome are "functional". Here we investigate the relative strengths of associations between coding and non-coding gene functionality and genomic features, by comparing reliably annotated functional genes to non-genic regions of the genome. We find that the strongest and most consistent association between functional genes and genomic features are transcriptional activity and evolutionary conservation. We also evaluated sequence-based statistics, genomic repeats, epigenetic and population variation data. Other features strongly associated with function include histone marks, chromatin accessibility, genomic copy-number, and sequence alignment statistics such as coding potential and covariation. We also identify potential issues with SNP annotations in short non-coding RNAs, as some highly conserved ncRNAs have significantly higher than expected SNP densities. Our results demonstrate the importance of evolutionary conservation and transcription activity for indicating protein-coding and non-coding gene function. Both should be taken into consideration when differentiating between functional sequences and biological or experimental noise.

Osteoporosis affects millions of women globally. In this study, we applied bioinformatics methods to screen for novel diagnostic biomarkers of osteoporosis in women using the GSE62402 and GSE56814 datasets. PCSK5, ZNF225, and H1FX were used to construct a diagnostic model. ROC, calibration, and decision curve analyses were performed to assess the diagnostic performance on the training (GSE56814) and external (GSE56815) datasets. The expression level of model genes was validated in GEO datasets. Furthermore, five transcription factors (ETS1, NOTCH1, MAZ, ERG, and FLI1) were identified as common upstream regulators of model genes. PCSK5, ZNF225, and H1FX serve as novel diagnostic biomarkers, providing new insights into the pathogenesis of and treatment strategies for osteoporosis in women.

The cockroach Blattella germanica possesses panoistic ovaries, in which oocytes lack nurse cells and therefore need to rely on their own transcriptional activity to support embryogenesis. Ovarian development in this species involves the development of a single basal ovarian follicle (BOF) per gonadotropic cycle, a process strictly regulated by endocrine signals, primarily juvenile hormone and ecdysone, which act at both the transcriptional and translational levels. In addition, transcriptional activity in these ovaries is necessary for both regulating and genome protection, and at this level, PIWI-interacting RNAs (piRNAs) play an essential role. Although insect ovaries are known to be particularly rich in piRNAs, their function in ovary maturation is still not well defined. For this purpose, we characterize the piRNA expression dynamics across seven key developmental and reproductive stages, ranging from late nymphal instars to post-vitellogenic adults. piRNA expression in B. germanica shows coordinated fluctuations. Expression remains stable in previtellogenic ovaries, whereas vitellogenic ovaries show pronounced changes. Moreover, vitellogenic ovaries exhibit reduced piRNA diversity due to strong enrichment of a subset of highly expressed piRNAs. Our data show that although piRNAs predominantly map to transposable elements, particularly LINEs, there is a notable increase in gene-derived piRNAs toward the end of the cycle. Our results suggest regulatory roles of piRNAs in modulating both TEs and mRNAs during BOF maturation, likely related to changes in the follicular cell program.

Yeasts ability to invade surfaces has important implications for infections and food contamination. Invasive growth in yeast is influenced by genetic and environmental factors. In this exploratory study, we investigated the effects of sodium sulfide, gene deletions, and environmental conditions on the invasive behaviour of the wine yeast strain AWRI 796. Sodium sulfide enhanced invasion in the (parent) AWRI 796 strain under nitrogen-limiting conditions, although its effect was obscured by experimental variability and pre-culture conditions. Genetic factors had a major effect on the overall invasive phenotype, with deletion of key genes suppressing invasion. Most gene-deletion mutants did not significantly affect how the colony responded to sulfide. In addition to sulfide and genotype, environmental conditions also influenced invasive behaviour. The pre-2xSLAD pre-culture condition was best for detecting sulfide-induced growth, and later plate washing time and decreased nutrient levels enhanced invasiveness. Our experimental design and findings provide a framework for understanding the determinants of yeast invasiveness, which may inform future studies on filamentous yeast behaviour.

Epigenetics may play a crucial role in livestock adaptation to environmental challenges like heat stress. In recent years, a growing number of studies have investigated the epigenetic mechanisms underlying dairy cow adaptation to heat stress. However, there is still limited knowledge about the effects of heat stress on immune cells and immune-related phenotypes. Herein we aim to identify heat-stress induced DNA methylation variations on blood methylome potentially affecting regulatory regions and associated phenotypes. Blood samples were collected and peripheral blood mononuclear cell (PBMC) isolated from four cows before (D0) and after (D14) a 14-d heat stress challenge (cyclical THI 72-82) and, from four cows kept in thermoneutral conditions (THI 61-64). Heat-stressed cows had ad libitum access to diets supplemented with adequate levels of vitamin D and Ca (12,000 IU/kg of vitamin D and 0.73% Ca, respectively). To eliminate confounding effects due to differences in nutrient intake, cows maintained under thermoneutral conditions were pair-fed (PF) to their heat-stressed counterparts and received adequate concentrations of vitamin D and Ca as well. Reduced representation bisulphite sequencing (RRBS) was used to profile PBMCs methylome. Differential methylation analysis was performed using methylKit and DSS softwares ({Delta}meth [&ge;] 25%, adjusted p-value < 0.01), retaining only commonly detected differentially methylated cytosines (DMCs). A total of 2,908 DMCs were identified when comparing pre- and post-heat stress samples. After excluding 649 DMCs that were also detected under thermoneutral conditions, as these changes were likely associated with feed restriction inherent to the pair-feeding design rather than with heat stress per se, 2,259 heat stress-specific DMCs remained, predominantly hypomethylated. About half of the DMCs are annotated in intronic and intergenic regions; known to harbor regulatory elements. By intersecting the DMRs with publicly available functional annotation data, we observed hypomethylation on regulatory regions putatively affecting cows immune system. As an example, we identified a loss of methylation within an enhancer region of the MSN gene, which is involved in lymphocyte homeostasis, and a loss of methylation in the promoter region of MECP2, a well-established epigenetic regulator with a central role in chromatin organization and gene expression. These findings highlight the impact of heat stress on dairy cow immunity and provide new insights into its epigenetic regulation under environmental stress. Interpretative summaryThis study examined DNA methylation changes induced by heat stress in dairy cows to elucidate epigenetic mechanisms of thermal adaptation. Using RRBS on PBMCs, 2,259 heat stress-specific differentially methylated cytosines were identified, predominantly hypomethylated and enriched in regulatory regions. Functional annotation highlighted immune-related pathways, including hypomethylated regulatory regions near genes (e.g., MSN, ZBTB33, SLC25A5, GNAS, FAM3A, and MECP2) associated with immune function. These findings indicate that heat stress induces targeted epigenetic modifications potentially affecting immune regulation in dairy cows.

BackgroundMaize knobs are regions of constitutive heterochromatin that are readily identified in both meiotic and somatic chromosomes. These structures have been characterized as stable throughout the cell cycle, exhibiting late replication during the S-phase, and are composed of two specific families of highly repetitive DNA sequences: K180 and TR-1. Although widely used as cytogenetic markers due to their variability in number and chromosomal position across inbred lines, hybrids, and landraces, little is known about their chromatin structure and dynamics. In this study, we analyzed chromatin accessibility of knobs using DNS-seq data across four maize tissues representing distinct developmental stages. ResultsOur results reveal that K180 knobs exhibit tissue-specific variation in chromatin accessibility, transitioning between open and closed states during development. In contrast, the TR-1 knob of chromosome 4 remained consistently inaccessible across all tissues analyzed. A knob composed of both K180, and TR-1 further supported this observation, with only the K180 region showing dynamic accessibility. To validate these findings, we also analyzed other repetitive regions such as centromeres, which showed a uniformly closed chromatin structure similar to TR-1. These results suggest a unique developmental modulation of chromatin accessibility associated with K180 repeats. While the chromatin accessibility of knobs does not reach the levels observed at Transcription Start Sites (TSS), the comparison among different classes of repetitive DNA within maize constitutive heterochromatin provides compelling evidence for sequence-specific and tissue-specific chromatin dynamics. ConclusionsOur findings uncover a previously unrecognized property of maize knobs and establish a reference for future studies on chromatin organization and epigenetic regulation of repetitive DNA in plant genomes.

Primarily inhabiting the harsh, high-altitude environment of the Qiangtang National Nature Reserve exceeding 5,000 meters above the sea (m.a.s.l.), the golden wild yak is critically endangered, with fewer than 300 individuals remaining in the world, a situation exacerbated by the significant challenges of conducting research and conservation of their genetic resources. Somatic cell nuclear transfer (SCNT) can be an effective method for their preservation, but facing several obstacles in this context, including the hypoxic stress at high altitude that impairs embryonic development due to in vitro manipulation, and constraints of long-distance embryo transport. In the present study, the ear tissue was collected from a childhood male golden wild yak at Xizang Geye Wildlife Rescue Station (4800 m.a.s.l.) and send to Institute of Animal Science at Beijing to derive fibroblast cells. Using fibroblast cells of the golden wild yak as nuclear donors, and bovine oocytes from a local slaughterhouse at Beijing as recipients, the interspecific SCNT (iSCNT) embryos were generated and in vitro developed to blastocysts. To maintain the embryonic viability after long-distance transportation from Beijing to Xizang, iSCNT blastocysts were subjected to cryopreservation by vitrification method. Thawing of vitrified iSCNT blastocysts were completed at Xizang Dangxiong Yak Breeding Innovation Base (4200 m.a.s.l.), and transferred into the uterine horn of domestic yaks. 257 days after blastocyst transfer, a cloned golden wild yak was successfully harvested on January 10, 2026. This work demonstrates, for the first time, that interspecies somatic cell nuclear transfer can successfully produce a cloned offspring under extreme conditions, spanning 4800 m.a.s.l. donor origin, long-distance vitrified embryo transportation, and high-altitude blastocyst transfer at 4200 m.a.s.l., establishing a viable strategy for conserving critically endangered high-altitude species.

BackgroundInterferon-gamma (IFN-{gamma}) is the primary effector cytokine of adaptive anti-tumor immunity, yet it paradoxically induces a potent immunosuppressive tumor microenvironment (TME). The full mechanistic scope of this paradox in head and neck squamous cell carcinoma (HNSC) has not been characterized at the transcriptomic scale. MethodsUsing TCGA HNSC RNA-seq data (n = 522), we applied an integrated computational pipeline: Spearman correlation analysis, principal component analysis (PCA), UMAP, K-means clustering (k = 4), Random Forest regression, deep neural networks, permutation importance, JAK-STAT cascade mapping, and DNN-based transcriptome-wide mediation analysis across 57 IFN-{gamma} pathway and 78 immunosuppressive genes. ResultsIFN-{gamma} pathway activity was universally and positively correlated with six immunosuppressive axes, including checkpoints (CD274; LAG3; IDO1), Tregs, myeloid suppression, and tryptophan catabolism. K-means clustering identified four immunologically distinct tumor subgroups. DNN models predicted suppressive TME. Permutation importance identified IRF8 as the dominant mediator linking IFN-{gamma} signaling to immunosuppression. DNN mediation analysis identified PDCD1LG2 (PD-L2) as the strongest intermediary between IFNG and PD-L1 regulation, followed by JAK2 and GBP5. ConclusionsIFN-{gamma} orchestrates coordinated immunosuppression in HNSC through JAK-STAT-IRF8 signaling. PDCD1LG2 and JAK2 are actionable mediators of this paradox, supporting combination strategies co-targeting IFN-{gamma}-induced checkpoint induction and direct checkpoint blockade in HNSC immunotherapy. GRAPHICAL ABSTRACT

BackgroundRed macroalgae (Rhodophyta) are ecologically and economically important marine primary producers, yet genomic resources for most species remain scarce. Delisea pulchra, a temperate red alga known for its halogenated furanone-based chemical defenses, serves as a model for studying algal-microbe interactions, antifouling mechanisms, and disease dynamics. ResultsHere we present a high-quality genome assembly of this species. The nuclear genome comprises 134 Mbp across 271 contigs with an N50 of 1.47 Mbp and encodes 13,387 predicted protein-coding genes. Comparative genomics with other red algae revealed expansions in gene families involved in DNA methylation, and oxidative stress responses, including glutathione S-transferases and superoxide dismutases. Analysis of glycosyltransferases, sulfatases, and sulfurylases implicated in galactan biosynthesis suggests D. pulchra possesses a complex and potentially novel extracellular matrix. We also identified several vanadium haloperoxidases (vHPOs), heme-dependent haloperoxidases (hHPOs), and two type III polyketide synthase (PKS) genes unique to the D. pulchra, which together represent promising candidate genes for bromofuranone production. ConclusionThe D. pulchra genome provides a foundation for molecular investigations into defense, signaling, and host-microbe interactions. It has been deposited at the European Nucleotide Archive under accession number PRJEB101077. All datasets, annotations, and interactive tools for exploring the genome are also available through the Rhodoexplorer portal at https://rhodoexplorer.sb-roscoff.fr.

This study investigated the mutagenic effects of ethyl methane sulfonate (EMS) on the M{square} generation in cowpea (Vigna unguiculata (L.) Walp.) cultivar Wang Kae. A total of 275 M{square} seeds were treated with EMS concentrations of 20 mM, 40 mM, and 80 mM (75 seeds per treatment) by soaking for six hours, while 50 untreated seeds served as the control (0 mM). Phenological, yield-related and yield traits were recorded, and data were analysed using Jamovi 2.7.15 and JASP 0.95.4.0 through one-way ANOVA with post hoc contrast, principal component biplot, and cluster analyses. No optimal mutagenic concentration (LD50) was identified. Seed germination and seedling survival rates increased with increasing EMS concentration, ranging from 70.00% and 62.00% in the control (0 mM) to 89.33% and 74.67% at 80 mM, following the trend 0 mM < 20 mM < 40 mM < 80 mM. Significant differences (P < 0.05) were observed among treatments for all phenological traits, pod length, locule number, seed traits, and yield per plant. Yield was significantly higher (P = 0.047) at 20 mM (61.19 {+/-} 3.34 g) compared to the control. Contrast analysis identified genotypes B33 and D56 as the most productive mutants, with yields of 125.44 g and 111.85 g, respectively. Principal component analysis extracted eighteen components, with the first four cumulatively explaining 50.60% of total variation. Biplot analysis of PC1 and PC2 captured all phenological traits, key seed traits, and yield attributes, highlighting the superior performance of B33 and D56. Cluster analysis partitioned the 190 genotypes into six groups, with B33 and D56 constituting distinct clusters. EMS mutagenesis effectively induced heritable phenotypic variation, with putative superior genotypes identified for advancement to M{square} and evaluation in replicated multi-environment trials toward the development of farmer- and consumer-preferred cowpea varieties.

Paralytic Shellfish Toxins (PSTs) are produced by certain species of cyanobacteria and dinoflagellates. Part of the PST biosynthetic pathway has been elucidated in cyanobacteria, and the implication of some sxt genes has been confirmed by experimental studies. Contrary to cyanobacteria, knowledge about PST biosynthesis in dinoflagellates is more limited and generally restricted to comparative studies with the cyanobacterial pathway. To investigate the specificity of the PST pathway in dinoflagellates, 16 toxic and non-toxic A. minutum strains from a recombinant cross were compared, without prior assumption on genes or metabolites involved in PST synthesis, using an integrative approach combining untargeted metabolomic and transcriptomic data. Among the 60 most distinguishing transcripts between toxic and non-toxic strains, only 3 sxt genes were present, sxtA4, sxtG, and sxtI. In contrast, non-sxt homologs were detected as highly discriminant between these two phenotypes. More specifically, a phyH homolog may act as the analog of sxtS found in cyanobacteria. Moreover, we identified four putative synthetic PST intermediates. Among these, Int-C2, correlated with the toxic phenotype, whereas 3 others were detected in both toxic and non-toxic strains, suggesting that these strains may share some parts of the biosynthetic pathway. Finally, our results showed that PST biosynthesis in dinoflagellate results from the activity of sxt genes, acquired by horizontal gene transfer from cyanobacteria, as well as from other genes not acquired from cyanobacteria, such as phyH.

It has been known for decades that bacteriophages encode tRNA genes, but their function and the factors contributing to their acquisition and retention are unclear. Although tRNAs are found in a variety of phages infecting a variety of bacteria, many large-scale computational studies investigating tRNA acquisition and retention in phages are specific to Mycobacterium phages; however, these findings may not be representative of other phages or bacteria. This work uses a broader sampling of phages and hosts to investigate the relationships between codon usage bias, infection cycle, and tRNA gene numbers in phage genomes. We analyzed 154 phages infecting 7 host genera, including Gram-negative (Escherichia, Shigella, Salmonella) and Gram-positive (Bacillus, Lactobacillus, Staphylococcus, Mycobacterium) bacteria. Phages included temperate and virulent representatives, plus a range of tRNA numbers and morphologies. All phages and hosts were analyzed using four metrics: GC content, Effective Number of Codons, Relative Synonymous Codon Usage, and tRNA Adaptation Index. On a global scale, virulent phages with many tRNA genes show greater differences in codon usage and codon adaptation compared to their respective hosts. Gram-negative bacteria and their phages generally exhibit greater differences in codon usage compared to Gram-positive bacteria and their phages. Phages infecting Gram-negative hosts also tend to encode more tRNA genes. In nearly all genus-level comparisons, Mycobacterium phages were different from any other host and from global patterns. This suggests previous computational studies performed in Mycobacterium phages are likely not applicable on a global scale or to phages infecting other host genera. AUTHOR SUMMARYBacteriophages, or phages, are viruses infecting bacteria. They are abundant in all environments, yet how they interact with their bacterial hosts is still not well-understood. Like other viruses, phages must rely on the host translational components to replicate and form new phage particles; and similarly to other parasites, phages have genomes that differ significantly from their hosts in terms of composition. In this work, we explore the relationship between phage lifestyle, number of tRNA genes encoded, and genome differences from the host using a variety of phages and their associated hosts. Phages can be either virulent (do not integrate into the host genome) or temperate (capable of integrating into the host genome), with differences from the host genome more pronounced in virulent phages. There are many phages that also carry tRNA genes, and having higher numbers of tRNAs is associated with larger differences from the host genome. The findings here indicate that virulent phages carrying large numbers of tRNAs diverge the most from host genome composition.

Text abstractsLipid homeostasis is essential for organismal physiology, and its disruption contributes to metabolic disorders. Using an unbiased genetic modifier screen in Drosophila, we identified GAR1, a core component of the box H/ACA small nucleolar ribonucleoprotein complex, as a pivotal regulator of systemic lipid storage. We show that the H/ACA snoRNP complex is essential for maintaining lipid droplet morphology in adipose tissue and preventing ectopic fat accumulation. Moreover, null mutants of Gar1 or Dkc1 exhibit severe developmental defects, including reduced body size and larval lethality. RNA-seq analysis revealed that Gar1 dysfunction triggered widespread alternative splicing defects, specifically targeting key transcripts within the insulin signaling cascade, including chico, Pi3K92E, sgg, and Lip4. Furthermore, knockdown of Gar1 impaired insulin signaling, as evidenced by the reduced membrane localization of the tGPH fluorescence. Genetic epistasis further positions GAR1 upstream of the lin-28/foxo axis, as knocking down lin-28 or foxo fully rescues the lipometabolic defects in GAR1-deficient animals. These findings reveal a previously unrecognized link between the snoRNP machinery and metabolic process, establishing the box H/ACA complex as an important coordinator that integrates RNA processing with insulin-mediated nutrient sensing to ensure developmental and lipid homeostasis. Article summaryLipid metabolism is tightly controlled by multiple factors. To find new regulators, the authors performed a genetic screen and identified a small nucleolar protein GAR1 participate in fat storage and larval development. They demonstrated a critical role of box H/ACA snoRNP complex in modulating alternative splicing and balancing insulin cascade. Blocking two insulin-related genes reversed the lipid defects caused by Gar1 loss. These findings revealed the box H/ACA complex integrates RNA processing with insulin-mediated nutrient sensing to ensure developmental and lipid homeostasis, offering a perspective for understanding the metabolic regulation network.

Phenotypic plasticity is a key mechanism by which plants adjust their traits to environmental changes. These phenotypic adjustments are driven by plastic changes in gene expression regulated by gene regulatory networks. Drought, a major selective force in Mediterranean ecosystems, provides a powerful context to examine how genomic plasticity translates into phenotypic responses. Here, we used Dianthus inoxianus, a drought-tolerant Mediterranean carnation, in order to characterize the phenotypic and transcriptomic plasticity in response to drought stress combining ecophysiological measurements with RNA-seq, gene co-expression and gene regulatory network analyses. Most of the phenotypic traits exhibited low plasticity in response to drought, except water and osmotic potential. At transcriptome level, we identified 57 plastic genes, suggesting that drought tolerance in D. inoxianus relies predominantly on constitutive gene expression. These plastic genes were enriched in processes typically related to drought response, such as cell wall components and abscisic acid (ABA) signaling. Some plastic genes belonged to drought-responsive modules, while others were hubs in different modules acting as inter-modular connectors. Furthermore, the regulatory network revealed that these plastic genes were strongly regulated by multiple stress-responsive transcription factors, and that drought-associated modules were regulated through both ABA-dependent and ABA-independent pathways. In addition, we identified contrasting patterns of canalization and decanalization, with immune and post-transcriptional regulation remaining canalized under drought, whereas photosynthesis and amino acid metabolism became decanalized, potentially releasing cryptic genetic variation. Overall, our results emphasise that drought tolerance in D. inoxianus emerges from a strategy combining preadaptation with targeted plasticity in key molecular pathways.

High and low tides occur twice a day (every [~]12.4 hours), with the largest tidal ranges during spring tides around new and full moons (every [~]14.765 days). While these lunar cycles are known to influence many animal phenotypes, particularly the reproduction of coastal animals, the genetic basis of lunar-related rhythms remains unclear. Since phenotypic variation is a valuable resource for elucidating such mechanisms, we examined geographic variation in the lunar-regulated mass spawning of the grass puffer (Takifugu alboplumbeus) along the Japanese coast. We found that western populations spawn during the first half of the spring tides, whereas eastern populations spawn during the second half. Furthermore, although spawning typically occurs a few hours before high tide, this timing is restricted to a specific time window that is earlier in the western populations than in the eastern ones. Behavioral analysis of larvae also revealed a shorter free-running circadian period ({tau}) in the western population than in the eastern ones. As differences in {tau} affect individual variation in the timing of physiological functions and behaviors, we hypothesized that differences in {tau} could account for the different time windows and consequently the observed difference in spawning days. Population genomics analysis identified proline-rich transmembrane protein 1-like (prrt1l) as a candidate gene. Expression of prrt1l was observed in the circadian pacemaker suprachiasmatic nucleus, and triple CRISPR F0 knockout of prrt1l shortened the free-running period in larvae. These findings suggest a potential mechanism underlying the geographic variation in lunar-synchronized spawning behavior. HighlightsO_LIThe geographic variation exists in the lunar-regulated spawning of the grass puffer, with differences in spawning dates and times between western and eastern Japan. C_LIO_LIThe free-running period of western populations is shorter than that of eastern populations, which is consistent with their earlier spawning timing. C_LIO_LIPopulation genomics analysis identified prrt1l as a candidate gene harboring population-specific missense mutations, the knockout of which shortens the free-running period. C_LI